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Abstract 

A theoretical technique which allows to include contributions from non- 
orthogonality of the electron states in the leads connected to a tunneling 
junction is derived. The theory is applied to a single barrier tunneling 
structure and a simple expression for the time-dependent tunneling current 
is derived showing explicit dependence of the overlap. The overlap proves 
to be necessary for a better quantitative description of the tunneling cur- 
rent, and our theory reproduces experimental results substantially better 
compared to standard approaches. 

PACS numbers: 72.10.Bg, 72.15.-v, 73.23.-b, 73.40.Gk, 73.63.-bk 

Achievements in nano-materials science is expected to have importance in many sci- 
entific fields, including information technology, quantum computing and fuel cells. In 
particular, tunneling phenomena have been under focus recently, both in magnetic het- 
erostructures as well as for quantum dot systems. The purpose of this paper is to develop 
an improved description of this phenomenon for general tunnel junctions, with possible 
application to the aforementioned scientific questions. 

To focus the discussion, we mention that conductance measurements on extremely 
small metal-insulator-metal (MIM) junctions were carried out by Vullers et a/.0 showing a 
non-linear conductance as a function of the bias voltage for low temperatures. The same 



behaviour has been reported for MIM double junctionsa and Ti/TiO x tunneling barrier 
systemsihi. The non-linearity in the current- volt age (J — V) characteristics appears for 
source-drain bias voltages larger than the spacing of the quasi one dimensional subbands 
since different numbers of subbands become available for transport in the forward and 
reverse directions!. In the study by Simmons@ the current was found to depend non- 
linearly on the voltage, roughly as V + 7V 3 . 

Many theoretical studies of transport in nanostructures with tunneling barriers rely on 
the transfer HamiltonianiHi! which contains serious inconsistencies^. The principle of the 
transfer Hamiltonian is a division of the system into subsystems. This is motivated by the 
fact that the physical properties of the subsystems may be different and, hence, require 
different descriptions. Another motivation is that one is directly offered the possibility to 
generalize the approach to any number of tunneling barriers in the system. The transfer 
(tunneling) between the subsystems arises due to an overlap of the wave functions in 
the region of the barrier whereas the electron operators of the different subsystems are 
assumed to be anti-commuting. Qualitatively this may be motivated since the leakage 
of a wave function in one subsystem into the other is exponentially small. The J — V 
characteristics given in this picture also shows a non-linear structure for large bias voltages. 
Quantitatively, though, the assumption of anti-commuting operators creates serious errors 
in the calculations of the current. This becomes particularly evident in the equilibrium 
situation displayed in Table |, in which the four lowest states of a particle in a one 
dimensional hard walled box with a scattering potential are given. The energy levels 
are, as expected, reproduced within the non-orthogonal representation (NOR) with much 
higher accuracy than in the orthogonal representation (OR). Attempts that go beyond 
the transfer Hamiltonian have been made, e.g. by expanding the non-orthogonal states 
into a new Hilbert spaceB. However, the proven success and physical transparency of 
the transfer Hamiltonian approach makes it desirable to extend its applicability to more 
general situations where the overlap is large, without making use of perturbation theory. 
This can indeed be achieved, which we demonstrate in this paper. 

In order to overcome the inconsistencies with the transfer Hamiltonian formalism, we 



develop a theoretical approach for time-dependent transport through tunneling systems 
in which the overlap between the subsystems give explicit contribution to the current. 
Technically, we will express the properties of the original system in terms of the operators 
constructed of the wave functions of each subsystem. The resulting model structurally 
resembles the transfer Hamiltonian, although the physical interpretation is different. We 
have chosen the single barrier system simply to show the features of our approach. The 
main result of this paper is the Eqn. @ for the time-dependent tunneling current through 
a single barrier. This expression is applied to a MIM junction in order to analyze the 
effect of overlap on the current. To our knowledge there does not exist any derivation nor 
analysis of time-dependent transport in tunneling junctions where the non-orthogonality 
is not disregarded. 

Let us now proceed starting with the one particle Hamiltonian 



where V is any potential describing a system of two leads with an insulating layer in 
between. We introduce the two potentials V a , a = L,R, for the left (L) and the right 
(R) subsystem, respectivelyi'SEl. For instance, the left potential can be written as Vj, = 
V(x)9(— x + ol) + V(ai)0(x — a^), where 6{x) is the Heaviside function and ol is a 
turning point for the left subsystem. In each subsystem there are orthonormal eigenstates 
{4>k,£ka}k*ea from which the corresponding field operator ip Q (t,x) = J2kaea Cka(t)<f>k(x) 
is constructed. Here t is time and x = (r, a) is a vector of the spatial coordinate r and 
the spin a. Suppose that ip is the field operator of the system formed by the potential 
V. Then, this operator can be expanded in terms of ip a by the trivial identity if)(t,x) = 
J2 a V'a^j x ) + [' l P(t, x )~J2 a "4>a(t, %)]■ Following reference^! we project ip onto the subsystem 



k G a, interpreted as the annihilation of a particle in the state 4>k with spin projection 
a. Creation c\ a of a particle in the state 4>k is defined similarly. These projections are 



H = £ + V, 



a by 
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possible to use directly for a second quantized form of the Hamiltonian. However, such an 
expansion gives an inconvenient expression of the Hamiltonian with the overlap matrix 
appearing explicitly. Thus, in order to proceed further, we define the operators 

c k *{t) = Ek>Okk>Ck>*(t), 

C L(*) = Sfc' (Cfefe 1 ') c\., a (t), 
where k' runs over all states in L U R and O kk , is the element kk' of the inverse of the 
overlap matrix of the wave functions (pk> given by Okw = (0fc|<?W) = Oy k . By a 
limitation to the case of spin conservation we can omit the spin indices in the overlap 
integral. The expectation value of the Hamiltonian in these operators is 

H = J ^H^dx = H L + H R + H T , (2) 

where we have defined TL a = J ip^Hipadx and Ht = J2aa'U iPaHip a 'dx + H.c). Here, we 
have neglected all expectation values that contain ip — J2 a Furthermore, we note that 
from the identity V = V a + [V — V a ], a = L,R, we find that the Hamiltonian of the lead 
1~La = T,kaea £ kacl a c k a + T,kk'ea(4>k\(y -Va^cpk^cl^Ck^, a = L, R. The last term is a sum 
of terms proportional to the integral of (pl(pk' over (or, oo) or (—00, ai) when a = L or 
a = R, respectively, in which domains the wave functions are exponentially small. Thus, 
this term is negligible and we arrive at the appealing form of the Hamiltonian 

7~t — ^ £paCp a c pa + ^ e qa c q(T c q(T + '^ / (v pg(7 c pa c qa + H.c), (3) 

paeL qaeR PQ<? 

where v pqrj = (4> p \H\<f) q ) is the mixing matrix element. The structure of the Hamiltonian 
(0) very much resembles the usual transfer Hamiltonian. Nevertheless, the meaning of 
the electron operators c ka , Ck a is altered, now carrying information of the full system 
rather than just of its subsystem. This fact is legible from the anti-commutation relation 
{cka,cl'a} = ®kk'- Indeed, when O^y ~^ $kk' we recover the transfer Hamiltonian with 
the usual interpretation of the operators Ck a . In this sense we conclude that the Eq. (|3]) 
generalizes the conventional transfer Hamiltonian. 

The expression in Eqn. @ is derived for the system in equilibrium. It is straight for- 
ward applicable to the non-equilibrium case by letting e fc(T — > Sk a {t) and v pqa — > v pqa (t). 
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For definiteness we derive an expression for the current flowing through the barrier from 
the left to right. The tunneling current through the barrier separating the leads is ex- 
pressed as the rate of change of the number of particles on, say, the left side of the junction 
(N L (t)) = J2pa( n P <r{t)) , where (n pa (t)) = {c pa (t)c pa (t)) . The time development of (n pa ) is 
given by the Heisenberg equation of motion yielding the tunneling current for each spin 
projection a 



Ut) = 2e3]T (V; qa (t)(cl(t)c p(7 (t)) +v; qa (t)O m 1 (cl(t)c p(7 (t)) 

pq \ ' 

= -2eK J2 (v; qa (t)F< a (t, t) - v; qa {t)0^g< a (t, t)) , (4) 



with the coefficients V pqa = v vqa + 0^ q e qa describing the tunneling. In Eq. (f|) we have 
identified the correlation function (c qa c pa ) with the lesser Green function F<(t,t) = 
i(c qa (t)c pa (t)) . This propagator is calculated within the non-equilibrium technique of 
Kadanoff and Baym0 for the Green function F pq(T (t,t') = (—i)(Tc pa (t)c qa (t')). From the 
equation of motion for F pqa (t,t') we obtain 

F pq(T (t,t') = g p(J {t,t')0; q l + f 1 g prj {tM)V pq(T {ti)g q a{tut')&t u (5) 

J 

where g ka = F kka is the conduction electron Green function satisfying the equation 
(id/dt — £ka)gka(t,t') = 5(t — t'). The contour integration in Eq. is brought to 
real time integration by the Langreth analytical continuation rulesEl, thus 



/oo r 
VpgAti) <&,(Mi)s£(*i,0 + »£(Mi)C('i>0 
-oo L 

The lesser, retarded and advanced expressions of the conduction electron GF are 

g k l(t,t')=if a (e ka )e- i f>^ dt \ 
g r £(t,t') = TiQ(±t T efc * (tl)dl 



dti. 



respectively, where f a (x) is the Fermi-Dirac distribution function. Before we con- 
tinue the derivation we rewrite the electron operators in terms of current states, i.e. 
c\ a (t) = cj^exp [i/j, a (t)], and c ka {t) = c^exp [—ifi a (t)]. This will explicitly show the ap- 
plied voltage dependence V(t) of the current, since Hl{P) — fin(t) = eV(t). Replacing the 



summation over p and q in Eq. (|4]) by energy integration in terms of the density of states 
p a (s a ) and noting that ^[(V^. — v* qa )0~ q l g^} = 0, the time-dependent tunneling current 
becomes 

J a (t) = -2e3? J V^(t)p a (e L )p a (e R ) 

x I* ^(ti)[/M-/(^)]e- l ^ (ey( ' 2)+(£i ^ ))d ' 2 dt 1 de i d^. (6) 

J — oo 

The mixing and the overlap are here replaced by the functions ViRait) = V a (eL, £r, t) and 
Olr = 0~ l (e L ,£ R ), respectively, satisfying V a (e pa ,e qcn t) = V pqa (t) and 0~ x [e V(J , e qa ) = 
0~ q . The formula ([]) reproduces results based on the transfer Hamiltonian in the limit 
of orthogonal subsystems, i.e. when Q^y — > 5kk'- It is important to note the fact that 
the tunneling coefficient Vlr^ = v LRa + 0^ R e Ra in our formulation, explicitly depends on 
the energies of the electrons involved in the conduction process. 

When V(t) = V and a stationary current is established through the barrier the Eq. 
(|) reduces to 

7r r w 

J a = 2e— j_ w \V LRa \ 2 [f(e - eV) - f(e))de. (7) 

This expression is given by assuming a constant density of states p a {e a ) = 1/2W, where 
2W is the conduction band width, and slowly varying mixing and overlap so that their 
respective values can be taken at the chemical potential, which are reasonable conditions 
for MIM-j unctions. In order to compare our theory with a realistic example we show 
in Fig. m the experimental J — V characteristics from Ref.ll (solid-dotted line) together 
with that of Eqn. (|7|) in both the non-orthogonal (solid line) and orthogonal (dashed 
line) representations. We have also included the corresponding result given by Simmons 
formula (dotted line)0. Note that Simmons formula and the orthogonal representation 
correspond to the standard methods used to calculate transport. From the figure, it 
stands clear that inclusion of the overlap contributes significantly to the behaviour of 
the J — V characteristics and the quantitative agreement is remarkably improved. The 
increase in the agreement with the experiments lies not only in the low voltage regime 
but also in that the current rises rapidly at a certain threshold voltage, which influences 
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the time-dependent current. For a 6 % increase in the barrier width our calculation 
(bold dash-dotted line) agrees exactly with the experimental results for positive voltages. 
The remaining discrepancy from the experimental curve, e.g. the observed asymmetry, is 
believed to stem from the lack of electron interactions in our model, for example charging 
effects. Moreover, in the simple calculations presented here we have merely computed the 
wave functions <p^ L and (f) flR , normalized to a unit probability at their asymptotic 

distances from the barrier x — > — oo and x — > oo. For simplicity we have used a rectangular 
potential barrier leading. 

In conclusion, we have developed a simple and transparent theoretical approach for 
time-dependent tunneling current through nanostructures which has a far wider appli- 
cability compared to standard methods. The ability of dividing the system into several 
subsystems, which then can be treated individually, is preserved without loss of accuracy 
when the inclusion of the overlap of the subsystems is allowed and all attractive features of 
the transfer Hamiltonian approach can be kept. The non-orthogonality is reflected in the 
non-zero anti-commutation relations of the electron operators of different subsystems. A 
formula for the time-dependent tunneling current through a single barrier structure, Eqn. 
([]), has been derived, which shows the necessity of including the overlap for a substan- 
tially better quantitative agreement with experiments. We also note that the formalism 
simply generalizes to the case of a two, or multiple, barrier structure. In particular, the 
region between the barriers can be interacting, for example a quantum dot. Then, a 
generalization to any number of contact leads is straight forward. 
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TRF) and the Swedish foundation for strategic research (SSF) are acknowledged. 
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The experimental data of the barrier's width and height given by Haraichi et alE are 
calculated from (3{V + 7 V 3 ) , which is only valid for low and intermediate voltages. 
In the Fig [I], we have used the more precise Jo[(<f> — eV/2) exp (—2s^2(j) — eV) — (0 + 
eV/2) exp (— 2s\/2</> + eV)]. The two formulae for the current are found in Ref.@. 

L.D. Landau and E.M. Lifshitz Quantum Mechanics 3 ed. (Pergamon Press 1977). 
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TABLES 

TABLE I. The four lowest energy levels of a 37 nm long hard walled box with a 5.3 nm wide 
and 178 meV high scattering potential located in the middle of the box. The energies (meV) 
are computed exact, with the overlap matrix taken into account (NOR) and ignored (OR). 



exact NOR OR 

20.265 20.266 18.866 

27.781 27.862 27.342 

83.868 83.592 79.383 

111.088 113.793 107.176 
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FIG. 1. The J — V characteristics of a 1.46 nm wide and 1.85 eV high MIM junction (the 
height measured from the equilibrium chemical potential) [19]. The experimental results by 
Haraichi et al. [5] (solid-dotted) is compared with the computations within the NOR (solid), 
NOR with a 6 % increase of the width (dash-dotted), OR (dashed) and Simmons formula 
(dotted) [7]. The equilibrium chemical potential is 1.75 eV and the conduction band width is 



2W = 40 eV. 
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